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Expressions for energy and angular momentum changes of the hydrogen atom 
due to interaction with the electromagnetic field during the period of the 
electron motion in the Coulomb field are derived. It is shown that only 
, the energy change for the motion between two subsequent passings of the 

pericenter is related to the quasiclassical dipole matrix element for transitions 
3 ' between excited states. 
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^ \ 1 Introduction 

Q ■ Classical hydrogen atom in a monochromatic electromagnetic field is one of the simplest 
■ real nonlinear system which dynamics may be regular or chaotic [Tl, l2], depending on the 
. relative field strength and frequency. Even one-dimensional classical model of highly excited 
atom yields results sufficiently close to the experimental findings. For theoretical analysis 
approximate mapping equations of motion, rather than differential equations, are most con- 
venient [2-7]. Here a two-dimensional map (for the scaled energy and for relative phase of 
^ ■ the field) is generalized for two-dimensional hydrogen atom, i.e. we calculate energy and 
angular momentum changes of the atom interacting with the electromagnetic field. 

2 One-dimensional atom in monochromatic field 

d . 

O ■ The Hamiltonian of the hydrogen atom in a linearly polarized monochromatic electromag- 
netic field (in atomic units) is jH [7] 
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Here P is the generalized momentum, c is the light velocity, 

cF 

A = sin(a;t + ?9) (2) 

is the vector potential of the field, F, uj and d are the field strength amplitude, field fre- 
quency and phase, respectively. The change of the electron energy can be obtained from the 
Hamiltonian equations of motion [8] 

E = -r ■ F cos(wt + 7?) . (3) 

One can introduce the scaled energy Eg = EjiJ^I^ and the scaled field strength Fg = 
Fjuj'^l'^. However, it is convenient [2-7] to introduce the positive scaled energy e = —2Es 
and the relative field strength Fq = Fg/el, with Eq being the initial scaled energy. 
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Integration of eq. ([3]) over the period of time between two subsequent passages of the 
electron at the apocenter results in the change of the electron energy [3l H] 



where 



ej+i = Ej — 7iFoelh{ej+i) siwdj , 
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Here s = e~'^l'^ = to / {—2EY^'^ is the relative frequency of the field, i.e. the ratio of the field 
frequency to the electron Kepler orbital frequency and J's{s) is the derivative of the Anger 
function. Introducing a generating function Glsj+ij'&j) [HI [IH] one can calculate the phase 
change over the period 

= + 2TTeJ^(^ - T[Foelri{ej+i) cos^j , 



where 
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Figure 1. Trajectories {e,'d) for the map ^ and ^ with the parameter nF^e^ ~ 0.0035 and initial condi- 
tions -do = TT, £0 = 0.3 - 0.003i {i = 0, 1, 2, . . .). 

Equations (jlj) and (jH)) describe the changes of the energy and phase in time. This map 
greatly facilitates numerical investigation of dynamics and ionization process. In figures [T] 
and figure [2] results of the numerical analysis of map (jlj) - ([6]) are presented. 



3 Two-dimensional atom in monochromatic field 

For calculation of the energy change of the arbitrary orientated two-dimensional atom in the 
electromagnetic field according to equation (3) one should perform the transformation of the 
coordinates (figure |3j). The change of the angular momentum of the atom it follows from the 
Hamiltonian equations of motion 

M = -— = -rF sin(a + cp) cos(wt + ?9) . (8) 
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Figure 2. Ionization tlircshold field dependence on the relative frequency sq. Continuous curve represents 
results calculated with variation of the initial phase -do, while dotted curve is for do = 0. 

By analogy with the scaled energy one can introduce the scaled angular momentum 
fjL — 2Ms — 2Mu^/^. Moreover, it is convenient to introduce the parametric equations of 

motion for the Coulomb potential [6-8]. 

Integration of equations (3) and (8) for the half of the period of the electron motion, i.e. 
for transition from apocenter to pericenter yield the scaled energy and angular momentum 
changes in the linearly polarized electromagnetic field 
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Figure 3. Two-dimensional atom in the electromagnetic field. 
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Here Es(2;)is the Weber function, z = es and e is the eccentricity of the elhpsis. By 
analogy or from equations (9) and (10) choosing appropriate initial phases of the field one can 
calculate the energy and angular momentum changes for the electron motion from pericenter 
to apocenter as well as for the complete period with different initial conditions. 

So, for motion between two subsequent passages at the apocenter we have generalization 
of eq. (4) for the energy change 
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and the angular momentum change 



/^i+i — H + 
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STT 
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One can calculate the energy and angular momentum changes for the electron motion 
between two subsequent passages at the pericenter in the similar way. 
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3.1 Approximation for relatively high frequency s 

For relatively high frequency of the field, s ^ 1, the asymptotic form of the Anger function 
Js(se) and its derivative J'^(se) may be used [7]. Then eq. (11) may be written in the form 



Ej+i = Ej + vrFoeo S (e - 2) 
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Introducing the generating function we can obtain the iterative equation for the phase dj as 
a generalization of eq. (6) 
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Using mapping equations (fT3|) and ( IT^ we can represent the energy and phase dynamics 
(figure H]) and calculate numerically the relative threshold field strength for the ionization of 
the two-dimensional hydrogen atom (figure |5]). 



3.2 Limiting cases for energy and angular momentum changes 
3.2.1 Approximations for very extended orbits 

For very extended orbits, e — )■ 1, the expansions for the energy ( fTTl) and angular momentum 
(1T2|1 in powers of /3 = -\/l — <^ 1 are 



Sj+i = Ej + ^^^^ {- (1 + Js(s) sin 79,- cos(^ 
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Figure 4. Phase plane (£,t?) of two-dimensional map ((T5|) and ((T^ with parameter ttFoEq = 0.0035, initial 
conditions 1)0 — tt, Eq — 0.3 — 0.003i (i — 0, 1, 2, . . .) and orientation angle (a.) Lp = and (b) 
(p — 7r/6. 
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For e = 1 eq. (15) coincides with eq. (4) for the one-dimensional atom, while eq. (16) 
represents change of the angular momentum, resulting to transition to the elliptic states with 
nonzero angular momentum. 



3.2.2 Approximations for almost circular orbits 

For almost circular orbits the eccentricity is small, e — )■ 0. Expansion of expressions ([11 
and ( |T^ in powers of e may be obtained from equations (3) and (11). The results are 
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Figure 5. Ionization threshold field dependence on the relative frequency sq for two-dimensional atom with 
the eccentricity e = 0.9 and orientation angle ip = (curve 1) and (p = tt/6 (curve 2). 
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These expansions are not valid for s = 1, 2, 3, 
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4 Atom in the circularly polarized field 

The equation for the energy change of the atom in the circularly polarized field is 

Ek = —eF {vx cos{u!t + -d) ±Vy sin(a;t + -i?)) . 



(19) 



Here sign " + " or " — " corresponds to the same and opposite directions of the electron and 
field rotations, respectively. 



For the electron, moving between two subsequent pericenter the energy and angular 
momentum changes are 
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It should be noted, that expression in the curly brackets in eq. ( I20l) coincides with the 
expression in the quasiclassical radial dipole matrix element in the velocity representation 
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This correspondence, however, takes place only for interaction of the hydrogen atom 
with the circularly polarized microwave field and for integration of the equations of motion 
between two subsequent pericenters. In general, the energy [31 [7| and angular momentum 
changes depend on the integration interval. So, for motion of the electron between two 
subsequent apocenters, i.e., the most distant from the nucleus points, where the electron's 
energy change is minimal, the energy change is described by the expression similar to eq. 
( I20I) but instead of J_5(es) and J'_g(es) we have the Anger function and its derivative of the 
positive order, J<i(es) and J'^(es). This interval has been used in P, HI |5] for derivation of 
the Kepler map for the one-dimensional hydrogen atom. 



5 Conclusion 

Analytical expressions for the energy and angular momentum changes of the two-dimensional 
hydrogen atom in linearly and circularly polarized electromagnetic fields are derived. It 
should be noted that in general the expressions are rather complicated. The approximate 
expressions for limiting cases of the parameters are more convenient for analytical and nu- 
merical analysis of the dynamics. 

The derived expressions are suitable for the three-dimension hydrogen atom as well, and 
may be generalized for analysis of the chaotic motion (due to the Jupiter perturbations) of 
comets and asteroids in the Sun system. 
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